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Abstract 

We establish Extended Thermodynamics (ET) of rarefied polyatomic gases 
with six independent fields, i.e., the mass density, the velocity, the temperature 
and the dynamic pressure, without adopting the near-equilibrium approximation. 
The closure is accomplished by the Maximum Entropy Principle (MEP) adopting 
a distribution function that takes into account the internal degrees of freedom of a 
molecule. The distribution function is not necessarily near equilibrium. The result 
is in perfect agreement with the phenomenological ET theory. To my knowledge, 
this is the first example of molecular extended thermodynamics with a non-linear 
closure. The integrability condition of the moments requires that the dynamical 
pressure should be bounded from below and from above. In this domain the system 
is symmetric hyperbolic. Finally we verify the K-condition for this model and show 
the existence of global smooth solutions. 


1. Introduction 


Rational extended thermodynamics [l| (hereafter referred to as ET) is 
a thermodynamic theory that is applicable to nonequilibrium phenomena 
with steep gradients and rapid changes in space-time, which may be out of 
local equilibrium. It is expressed by a symmetric hyperbolic system of field 
equations with the convex entropy. 

As ET has been strictly related to the kinetic theory with the closure 
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method of the hierarchy of moment equations, the applicability range of the 
theory has been restricted within rarefied monatomic gases. Only recently, 
however, the ET theory of dense gases and of polyatomic rarefied gases 
has been successfully developed and has been obtained a 14-field theory 
that, in the limit of small relaxation times (parabolic limit), reduces to the 
Navier-Stokes-Fourier classical theory [^. This new approach to the case 
of polyatomic rarefied gases, in particular, is in perfect agreement with the 
kinetic theory in which the distribution function depends on an extra variable 
that takes into account the internal degrees of freedom of a molecule [^. 

Precise modeling of polyatomic gases and of dense gases in nonequilib¬ 
rium is an active and urgent issue nowadays with many important appli¬ 
cations like the study of shock wave structure , which is essentially 

important, for example, for the atmospheric reentry problem of a space ve¬ 
hicle. 

There are at least three different methods of closure of moment theory 
associated with the Boltzmann equation: 


1. A closure at the kinetic level proposed firstly by Grad in the case of 13 
moments, which is based on a perturbative procedure of the distribution 
function in terms of the Hermite polynomials Q; 


2. The phenomenological closure of ET by using the univers^ principles 
of physics to select admissible constitutive equations [l|, 


3. The kinetic closure of molecular ET by using the Maximum Entropy 
Principle (MEP) 


1C 


m. 


It is very suitable and in some sense surprising that the three different clo¬ 
sure methods give the same result in the case of the 13-moment theory for 
monatomic gases [l| and also in the 14-moment theory for rarefid polyatomic 
gases [^, provided that the thermodynamic processes are not far from equi¬ 
librium. 


We here want to focus mainly on MEP and therefore we want firstly to 
summarize the principle and its limitation. The principle of maximum en¬ 
tropy has its root in statistical mechanics and is developed by E. T. Jaynes 
in the context of the theory of information related to the concept of the 
Shannon entropy 13, ll^ . MEP states that the probability distribution that 
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represents the current state of knowledge in the best way is the one with 
the largest entropy. Concerning the applicability of MEP in nonequilibrium 
thermodynamics, this was originally motivated by the similarity between the 
field equations in ET and the moment equations, and later by the observa¬ 
tion made by Kogan 1^ that Grad’s distribution function maximizes the 
entropy. The MEP was proposed in ET for the first time by Dreyer [l^. The 
MEP procedure was then generalized by Muller and Ruggeri to the case of 
any number of moments 1]|, and later proposed again and popularized by 
Levermore 151]. In the case of moments associtaed to the Boltzmann equa¬ 


tion the complete equivalence between the entropy principle and the MEP 
was finally proved by Boillat and Rug geri 16l |. Later MEP was formulated 
also in a quantum-mechanical context 


17 


m 


As seen below, the truncated distribution function of MEP has the mean¬ 
ing also far from equilibrium provided that the integrals of the moments are 
convergent. The problem of the convergence of the moments is one of the 
main questions in the far-from-equilibrium case. In particular, the index of 
truncation of the moments N must be even. This implies that the Grad 
theory with 13 moments is not allowed in a situation far from equilibrium! 
For this reason the truncated distribution function is formally expanded in 
the neighborhood of equilibrium as a perturbation of the Maxwellian distri¬ 
bution. All closures by the MEP procedure are valid only near equilibrium. 
As a consequence, hyperbolicity exists only in some small domain of config¬ 
uration space near equilibrium [H, 19|. 

The aim of this paper is to prove that, in the case of rarehed polyatomic 
gases, a theory can be established with the closure that is valid even far from 
equilibrium. This is a theory with 6 independent fields, i.e., the mass density, 
the velocity, the temperature and the dynamic pressure. We will show that 
this non-linear closure matches compRtely the previous result obtained by 
using only the macroscopic method 


20 |]. 


2. Rarefied polyatomic Gas 


A crucial step in the development of the kinetic theory of rarefied poly¬ 
atomic gases was made by Borgnakke and Larsen 2l|. It is assumed that 
the distribution function depends on, in addition to the velocity of particles 
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c, a continuous variable I representing the energy of the internal modes of 
a molecule. This model was initially used for Monte Carlo simulations of 
polyatomic gases, and later it has been applied to the derivation of the gener¬ 
alized Boltzmann equation by Bourgat, Desvillettes, Le Tallec and Perthame 
22l |. The distribution function f(t, x, c, I) is defined on the extended domain 
[0, oo) X X X [0, oo). Its rate of change is determined by the Boltzmann 
equation which has the same form as in the case of monatomic gases: 


dtf + Ci dif = Q, 


( 2 . 1 ) 


where the right-hand side, the collision term, describes the effect of collisions 
between molecules. The collision term Q{f) now takes into account the 
existence of the internal degrees of freedom through the collisional cross 
section. Here dt = d/dt and di = djdxi. 


The idea, firstly proposed at the macroscopic level by Arima, Taniguchi, 
Ruggeri and Sugiyama ^ and successively in the kinetic framework in 
and 2^, is to consider, instead of the typical single hierarchy of moments, a 
double hierarchy, i.e., an T-series at the index of truncation N and a G-series 
at the index M: [N, M) system given by 


dtF + diFi = 0, 

(9tFfci + diFik^ = 0, 

dtFkik2 + diFikik2 = Pkik2-, dtGii + diGiii = 0 , 

: dtGiij^ + diGiiij^ = Quj^i 


dtFk]^k2...k]si T diFikik2...kff — Pkik2...kff) 


dtGn 


JIJ2---3M 


+ diGii 




= Qii 


J1J2---JM ■ 


with 


/* 1*00 

Fkik 2 -kp= / mf{t,x,c,I)ckiCk 2 ---Ck„‘p{I)dI dc, (2.2) 

Jr3 Jo 


Gllkik2-"kq 


[ [ mf{t,x,c,I)(c^+ 2 —]ck^Ck 2 ---Ckq(p{I)dIdc, 

JmJo V 


(2.3) 
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and 0 < p < N, 0 < q < M (when the index p = 0 we have F and 
when q = 0, Gu). The double hierarchy is composed of the traditional 
velocity-moments F’s and the energy-moments G’s where the variable I of 
the internal modes plays a role. The connection between the index M and 
N is discussed in 23|. The non-negative measure (p{I) dl is introduced so 


as to recover the classical caloric equation of state for polyatomic gases in 
equilibrium. The functional form of ip will be given in the next section. 


2.1 Equilibrium Distribution Function and the Euler System 

Let us consider firstly the case of 5 moments corresponding to an Euler 
fluid. In this case, = 1 and M = 0. The collision invariants in this model 
form a 5-vector: 

m(l,ci,c^ + 2-) , (2.4) 

\ my 

which leads to hydrodynamic variables: 




p \ 

r f°° 

( 

1 \ 

F. 

= 

pVi 

= J ^ 


Ci 

\Gu) 


pv"^ -L 2pe j 

Jr? Jo 


F21jm j 


f{t,x,c,I)ip{I)dI dc, 
(2.5) 

The symbols are the usual ones: p,Vi,e are, respectively, the mass density, 
the i-th component of the velocity and the specific internal energy. The 
entropy is defined by the relation: 


'R3 Jq 

By introducing the peculiar velocity: 


/• roo 

h = -kB / f log fip{l)dl dc. 
Jr3 Jo 


( 2 . 6 ) 


Gi — Ci Vi , 


(2.7) 


we rewrite Eq. (12.5p as follows: 



5 









Tommaso Ruggeri 


Note that the internal energy density can be divided into the translational 
part peT and the part of the internal degrees of freedom psf. 

per = f [ ^mC‘^f{t,:si,C,I)(p{I)dIdC, 

JR^ Jo ^ 

p poo 

pei= / Ifit,^,C,lMl)dIdC. (2.9) 

Jr? Jo 

The energy pex is related to the kinetic temperature T\ 

er = |—T, (2.10) 

2 m 

where ks and m are the Boltzmann constant and the atomic mass, respec¬ 
tively. The weighting function p{I) is determined in such a way that it 
recovers the caloric equation of state for polyatomic gases. If D is the de¬ 
grees of freedom of a molecule, it can be shown that the relation p{I) = 
leads to the appropriate caloric equation of state: 


e = 


O ’ 

2 m 


a = 


D-5 

2 


( 2 . 11 ) 


The maximum entropy principle is expressed in terms of the follow¬ 
ing variational problem: determine the distribution function /(t,x, C,/) 
such that h max, under the constraints (12.5p . or equivalently, due to the 
Galilean invariance, under the constraints (12.8p . The result due to Pavic, 
Ruggeri and Simic is summarized as follows: 

Theorem 1. The distribution function that maximizes the entropy 
under the constraints 112. <§)) has the form: 

“ m{kBTy+^T{l + a) ^ ‘ 

( 2 . 12 ) 


This is the generalized Maxwell distribution function for polyatomic 
gases. In [^, the following theorem was also proved: 

Theorem 2. If 112.12\) is the local equilibrium distribution function with 
p = p{t,x), V = v(t,x) and T = T(t,x), then the hydrodynamic variables p, 
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V and T satisfy the Euler system: 

dp d . . 


d d 


(2.13) 


+ pe + py+p|»4=0 


with 


rj. 

P = — pT, 
m 


D ks 

e = - -T. 

2 m 


(2.14) 


This is an important result because we can obtain the Euler equa¬ 
tions from the kinetic equation for any kind of polyatomic gases as well 
as monatomic gases. 

3. The 6 Moment-Equations for Polyatomic Gases 


The 14-field theory, N = 2 and M = 1, gives us a complete phenomeno¬ 
logical model but its differential system is rather complex and the closure 
is in any way limited to a theory near equilibrium. Let us consider now a 
simplified theory with 6 fields (referred to as the ET6 theory): the mass 
density p, the velocity v, the temperature T, and the dynamic (nonequi¬ 
librium) pressure H. This simplified theory preserves the main physical 
properties of the more complex theory of 14 variables, in particular, when 
the bulk viscosity plays more important role than the shear viscosity and 
the heat conductivity do. This situation is observed in many gases such 
as rarefied hydrogen gases and carbon dioxide gases at some temperature 
25, [^. ET6 has another advantage to offer us a more affordable 


ranges 


24 


hyperbolic partial differential system. In fact, it is the simplest system that 
takes into account a dissipation mechanism after the Euler system of perfect 
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fluids. In the present case we have: 
dF dFk 


dt dxk 

OFi ^ dFik 


dt dxk 
dFg dFuk 
dt dxu 


= 0 , 

= 0 , 

= Pii, 


with 


/ f\ 


( " ^ 

F^ 

= 

PVi 

\Fg) 


y pv"^ + 3{p + n) y 


dGii dGiik _ 

dt dxu 


I f 

Jr3 Jo 


( 1 


m 


a \ fIodide 


and 


/* 1*00 

Gg = pv‘^ + 2pe = / m(c^ + 2//m)//“ dl dc. 

Jr3 Jo 


(3.1) 


(3.2) 


(3.3) 


3.1 Nonequilibrium distribution function 

We want to prove the following theorem: 

Theorem 3. The distribution function that maximizes the entropy Ii2.0\) 
under the constraints 113.^} 113. 31) has the form: 


3/2 


1 + 0 


/ = 


m 


m {kBTy~^°‘T{l + a) \^2'KkBT 1 + ^ 


1 - 


3 n 

2(l+o) p , 


1 


exp 


\ 2 


\mC^ 


1 + 


n 


+ I 


1 - 


3 n 

2(l+a) p , 


(3.4) 


All the moments are convergent provided that 


n 2 

—1 < — < -(1 + a), a > —1. 
p 3 


(3.5) 


Proof: The proof of the theorem is accomplished with the use of 
the Lagrange multiplier method. Introducing the vector of the multipliers 
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(A, Xi, Xu, fill), we define the functional: 


C = - 


/ /' 

Jr3 Jo 


ksf log / /" dl dc -L A /9 - 


/_R3 Jo 


' \ 
mf Iodide] + 


+Xi (^pvi - L[ mfci /" dl dc^ + Xu ^pv'^ + 3(p + 11) - If mc^f - 




^pv'^ + 2p£ 



As this is a functional of the distribution fnnetion / and we want to maximize 
it with respect to / with the given macroscopic quantities, this functional 
can be substituted by the following one: 


/* roo r roo 

C = — / fes/log//" d/dc — A / / mf 1°" dl dc— 

Jr? Jo Jr? Jo 

p poo p poo 

—Xi / / mfci /" dl dc — Xu / mc^f /" dl dc — 

Jr^ Jo Jr? Jo 

-Wi/ay m + 2^ f dl dc. 


(3.6) 


m 


Since £ is a scalar, it must retain the same value in the case of zero hydro- 
dynamic velocity v = 0 due to the Galilean invariance. Therefore: 


C 



A:B/log/rd/dC 


p poo 

X / m//"d/dC- 

Jr3 Jo 




fl^ dIdC. 


mC^fl^dldC- 


(3.7) 


Comparison between (j3.6h and (13.7h yields the relations between the La¬ 
grange multipliers and the corresponding zero-velocity Lagrange multipliers 
indicated by hat: 


A — A — XiVi -|- {Xu -L pu)v‘^', Xi — Xi — 2{Xu + fiu)vi', Xu — Xu pu — Piu 

(3.8) 

which dictate the velocity dependence of the Lagrange multipliers. We no¬ 
tice that these relations are in accordance with the general results of the 
Galilean invariance 2^. The Euler-Lagrange equation dCjSf = 0 leads to 


"d/dc 
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the following form of the distribution function: 


= exp ''s , 


(3.9) 


where 


X — A + XiCi + XiiC'^ + fill ( C‘^ + 2 


I 


m 


By introducing the following variables: 


TTl '' Tfi 2 

^ = -;—{Xii + flu), Vi = 1 —Aj, C = 1 —Wi, ^2 = exp 

KB Kb kb 


-1-^X 

Kb 


(3.10) 


the distribution function can be rewritten as 




(3.11) 


Inserting (I3.11|) into the second equation of (j3.2p evaluated at the zero ve¬ 
locity, we obtain immediately = 0. Then the remaining equations of (13.2p 


and ()3.3p evaluated for v = 0 become 

r roo 

p= / mfl°' dl dC = m7r^/^r(l + a) 

Jr3 Jo 


n 


^3/2^1-hc 


1 P 

/ f mf(ff + —)l‘‘dIdC^ 

Jr3 Jo 


P + n = 

pe = 

Ir3 Jo 

= m7r^/^r(l + a) 


2 m 

4^5/2^1+a (^ + ^ C J ■ 


From the integrability condition, we have 


C > 0, > 0, a > -1. 


(3.13) 
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From (j3.12p and (j2.1ip . we obtain 




m 


m ' C j ’ 
3/2r f ^ “ 3\ Q 


\ 2 J ^5/2(A 




n = 


TTh Q /9 ^ J 

— TT^/^r - 


D-3\ D-3 1- ^ 


15 


m ^ 


2 \ 2 J D ^5/2^ 

We can invert these relations as follows: 


D-3 

2 


Q. 




P 1 

2pl + f ’ 


^ ^ ^ (I) - 3) 


m 2/95 — 3(p -|- n) mp 1 — 


= 


P 1 

n 

p 


D-3 

2 


m7r3/2r (-^^) l2pl-LSy Impl- 


3 n 


D-3 p 


(3.14) 


(3.15) 


The integrability conditions (j3.13p imply that, for a bounded solution, 
the ratio H/p must satisfy 


p -|- n > 0 and {D — 3)p — 311 > 0. 


(3.16) 


Inserting (j3.15p into the distribution function (j3.1ip . we obtain (13.4h 
and the proof is completed. When IT ^ 0 the (j3.4p becomes the equilibrium 
distribution function ()2.12p . 


3.2 Closure and Field equations 

Substituting (13.4p into the fluxes Fuk, Guk and into the production term 
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Pii of (13.111 . we obtain after some calculations 

p poo 

Fik= / mciCkfl'^ dl dc = pviVk + {p + 11)5**., 

Jq 

p poo 

Fiik= / / mc^ckfl°'dldc= [b{p + Yi)+ pv^)vk, 

Jb? Jo 

Guk = J ^ j Cfc//" dl dc = {pv^ + 2pe + 2p + 2U)vk, 

p poo 

/ / mC^Qif). 

Jr3 Jo 


Pii = Pii = 


i/“ dIdC. 


(3.17) 


From the balance equations of momentum and of energy in continuum 
mechanics, we know that 

Fik = pviVk - tik, Giik = {pv^ + 2pe)vk - 2tikVi + 2qk, 


where tik = —pdik+(^ik is the stress tensor, fT**, = —n5jfc+(T<*fc> is the viscous 
stress tensor, 11 is the dynamical pressure (non equilibrium pressure), iT<*fc> 
is the shear stress tensor that is a deviatoric tensor (traceless) and qk is 
the heat flux. Comparing with the closure p.l7ll (x,3)- conclude that the 
closure gives a result that in the 6-moment theory a<^ik> = 0 and qk = 0. 
This is the expected result that there exist no shear viscosity and no heat 
conductivity in the 6-moment theory. For what concerns ()3.17p 9. taking into 
account the Galilean invariance with c* = Gi + Vi, we obtain the zero-velocity 


of Fiik: 


Fiik = 


7r3 Jq 


mC^Ckfl'^dldC = 0 . 


Concerning the production term (j3.17P /i . the main problem is that, in order 
to have explicit expression of the production, we need a model for the colli¬ 
sion term, which is, in general, not easy to obtain in the case of polyatomic 
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gases. With (j3.17p we obtain the differential system of 6 moments: 


dp d,. 

Tt + = "• 

^ Kr + = “• 

d d 

— (2pe -L pv'^) + ^ { [2(p + n) -L 2p£ + pv"^] Vi] = 0, 

[3{p -L n) -L pv"^] + ^ { [5(p + n) -L 


(3.18) 


dt 


This system with the thermal and caloric equations of state (I2.14p is a closed 
system for the 6 unknowns {p,Vi,T,Il), provided that we know the collision 
term in (l3.17P /i . These results are in perfect agreement with the results 
derived from the phenomenological theory 20|, where a possible expression 
of the production term Pa is adopted. In the case of the BGK approach we 
obtain Pu = — 3n/T. 


3.3 Entropy density 


Let us study the entropy density h with non-linear distribution function: 

h = -kBj y^/log//"d/dC = - logll^ , 

with If given by (j3.15p q. The equilibrium part of the entropy density /ig is 
expressed as 

-log^Ey (3.19) 

Moreover we may notice that the chemical potential g = £ + ^—Ts is 
expressed as 

= :^(l + logL!s). (3.20) 

I m 

On the other hand, the non-equilibrium part of the entropy is expressed 


ks 


He = ps = —p — 


m 


D 
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as 


k = -{h - He) = log 
p m He 


(3.21) 


Since 


n 

VIe 


= 1 + 


k is expressed as 


, ks 


n 

1 H- 

P 


1 - 


D-3 

3 nN 2 


D — 3 p 


(3.22) 


1 - 


3 n 


D-3 


(3.23) 


D — 3 p 

This expression also coincides with the one obtained by the phenomenolog¬ 


ical ET approach 20|. k depends on a single variable Z = II/p; for D > 3, 

k exists and is bounded in the domain that contains the equilibrium state: 

D-3 


—p<Il<p- 


(3.24) 


in which k{Z) < 0,VZ ^ 0 and k has a global maximum A;(0) = 0 at 
the equilibrium state. Therefore the convexity condition is satisfied and the 
ratio between the dynamical pressure and the equilibrium one II/p satisfy the 


inequalities (|3.24p . According to the Theorem proved in 271], the entropy 
h has, as is expected, the maximum value at the equilibrium state where 
h = ps. 


3.4 Main Field and Symmetric System 

Taking into account (|3.15p , (|3.inp and (|3.8h , we obtain the full expression 
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of the Lagrange multipliers after some cumbersome calculations: 


g ks n / n 

-log —-h — I 1 H- 

T m Ue 2T \ p 


P 

Pll = — { I - 


-1 


2T 


D — 3 p 


(3.25) 


1 D u f n 

“ “ D-3'p V P 


-1 


1 - 


3 n 
D — 3 p 


-1 


According to the general theory developed in 16( for monatomic gases 
and in 23( for polyatomic gases, the components of the Lagrange multipliers 
coincide with the components of the main field 28| for which the original 
system ()3.18p becomes symmetric hyperbolic. Notice that, in equilibrium 
where IT = 0, the first five components of the main field (j3.25p coincide with 
those obtained by Godunov for the Euler fluid 29|: 


Me = —; 


M\e = 


PlllE = 


2r’ 


while Xii\e = 0 according to the fact that the Euler fluid is a principal 
subsystem of the 6-moment system 30(. In the reference 20|, we proved 
that the system can be equivalent to the one obtained many years 


ago by Meixner 
paper 


31 


32l ] via the internal-variable procedure. In the same 
20|, it was also proved that, in the limit of monatomic gas D —3, 


the system (j3.18p converges to the Euler system provided that initial data 
are compatible with the case of monatomic gases, i.e., we should choose 
n(0, x) = 0. 


We observe that, if we apply the so-called Maxwellian iteration 33|, the 
last equation of (j3.18p with the production given by the BGK approximation 
reduces to the Navier-Stokes equation in the absence of the shear stress 

3, 


35 


20 (]: 


n = —z^divv, with v = 


2D-3 

E =-7-—pr, 

3 L» 


(3.26) 


where v is the bulk viscosity. The system (|3.18l) in which the last equation is 
replaced by (I3.26h was studied by Secchi [361] and by Erid and Shelukhin [33, 
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38l |. Shock wave structure was recently studied for this model by Taniguchi, 
Arima, Ruggeri and Sugiyama showing an excellent agreement with ex¬ 
perimental data. 


3.5 K-condition, acceleration waves, and global smooth solutions 


We now want to prove that the so-called Shizuta-Kawashima K-condition 
39[] is satisfied by our differential system and therefore, according to the 
general theorems [40|, |4l|, |27t], in contrast to Euler fluid, there exist global 
smooth solutions provided that initial data are sufficiently smooth. The 
system (j3.18l) is a particular case of a generic system of balance laws: 


dFO(u) (9F*(u) 
dt 


(3.27) 


We recall that the system (j3.27p satishes the K-condition if, in the equi¬ 
librium manifold, any right characteristic eigenvectors d of (I3.27h are not in 
the null space of Vf, where V = d/du: 


Vf d|^ ^ 0 Vd. 


(3.28) 


Lou and Ruggeri 4^ noticed a connection between the K-condition and the 
global existence of acceleration waves and they propose a necessary weaker 
K-condition requiring (|3.28l) only for the right eigenvectors corresponding to 
genuine non linear waves. 


For a quasi-linear hyperbolic system, it is possible to consider a par¬ 
ticular class of solutions that characterizes the so-called weak discontinuity 
waves or, in the language of continuum mechanics, acceleration waves. Let 
us study a moving surface (wave front) T prescribed by the Cartesian equa¬ 
tion t) = 0 that separates the space into two subspaces. Ahead of the 
wave front we have a known unperturbed field uo(x,t), and behind an un¬ 
known perturbed field u{x,t). Both the fields uq and u are supposed to be 
regular solutions of (j3.27j) and to be continuous across the surface T, but to 
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be discontinuous in the normal derivative, i.e., 


[u] = 0, 


5u 

dff) 


= n^o, 


(3.29) 


where the square brackets indicate the jump at the wave front. In 4^ it was 
verified that the K-condition is equivalent to the relation: 


5fb = (Vf • H b oc (Vf • d) b / 0, 


where the operator <5 is dehned by <5 = \d/d(j)]. By introducing the mate¬ 
rial time derivative, the system (I3.18p in the BGK approximation can be 
rewritten as 


dvh 

P + P-^ = 

M + ^(p + n) = 0, 

p£ + {p + Ii)^ = 0 , 

OXk 


/ p + n 2 n 

\ P 3 y pr 


(3.30) 


As is well known, the characteristic velocities U and the right eigenvectors 
can be obtained from the system (13.271) by utilizing the chain rule of the 
operators: 

d d 

— —U6, — - Uid, f —>■ 0, 

Ot OXi 

and, in particular. 


• V = U-Vn, Vn = Vini. 

In the present case, from the system (I3.30p . we obtain 

-V5p + p5vn = 0 , 

— pFdv -L n(I(p -|- n) = 0, 

- pVbe + {p + n)(Iun = 0, 


- V5 


^-^61=0. 
P 3 


(3.31) 
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Taking into account the constitutive equations (I2.14jl and evaluating (j3.81ji 
in an equilibrium state, we have 


1) V = 0 <—> U = Vn, contact waves, (3.32) 


with Sp, Svt, Sp arbitrary (multiplicity 4) and Svn = 0, 511 = —6p (v-p 
denotes the tangential velocity); 


2 ) 



with 6p arbitrary, 


U = Vn± 



sound waves. 


5v = nV 


Sp 


Se = —Sp, 

Dp 


6U = 


3T)2 


(D - 3)eSp. 


(3.33) 


We notice that the sound velocity in (j3.33l) is independent of the degree of 
freedom D and coincide with the sound velocity of monatomic gas. This 
curious fact was explained by a general theorem 23|] in which was proved 
that for particular choice of (N, M) systems in which belong the 6 moment 
theory the characteristic velocities are independent on D. 


As only the last component of the production term f of the generic 
system (|3.27p is non-zero (see (|3.30p i. the K-condition (|3.28p is satisfied if 
5n 7^ 0. This is true for contact wave and for sound waves. Therefore the 
K-condition is satisfied and, together with the convexity of the entropy, we 
can conclude that, according to the general theorems, the 6-moment system 
has global smooth solutions for all time and the solution converges to the 
equilibrium one provided that the initial data are sufficiently smooth. 


4. Conclusions 


In the present paper, we deduced the system of equations for a dissi¬ 
pative fluid in which the dissipation is due only to the dynamical pressure. 
The closure was obtained by the method of the Maximum Entropy Principle 
without assuming that the processes are near equilibrium. This system is 
the simplest example of non-linear dissipative fluid after the ideal case of 
Euler. The system is symmetric hyperbolic with the convex entropy den¬ 
sity and the K-condition is satisfied. Therefore, in contrast with the Euler 
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case, there exist global smooth solutions provided that the initial data are 
sufficiently smooth. The result obtained here is in perfect agreement with 
the one obtained by using only phenomenological theory of ET 20|. The 
comparison with experimental data in the case of shock waves is excellent 
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